Quantitative Theory of a Time-Correlation Function in a One-Component 
Glass-Forming Liquid with Anisotropic Potential 
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The Shintani-Tanaka model is a glass-forming system whose constituents interact via anisotropic 
potential depending on the angle of a unit vector carried by each particle. The decay of time- 
correlation functions of the unit vectors exhibits the characteristics of generic relaxation functions 
during glass transitions. In particular it exhibits a 'stretched exponential' form, with the stretch- 
ing index ft depending strongly on the temperature. We construct a quantitative theory of 
this correlation function by analyzing all the physical processes that contribute to it, separating 
a rotational from a translational decay channel. Interestingly, the separate decay function of each 
of these processes is temperature independent. Taken together with temperature-dependent 
weights determined a-priori by statistical mechanics one generates the observed correlation function 
in quantitative agreement with simulations at different temperatures. This underlines the danger of 
concluding anything about glassy relaxation functions without detailed physical scrutiny. 
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The glass-transition is often described as 'mysterious' 
Q]; it is the opinion of the present authors that there is 
nothing mysterious in the glass-transition, it is only nec- 
essary to penetrate the inhomogeneous states that de- 
velop naturally and understand the role of a variety of 
dynamical processes that appear in various examples of 
glass transitions. Different systems may exhibit differ- 
ent ways of glassifying, although there are many generic 
aspects to identify the phenomena as belonging to the 
same class. To underline this view we present in this 
Letter an in-detail analysis of one very interesting model 
of glass-formation, exposing how a number of distinct 
process come together to exhibit a seemingly complex 
phenomenology. In particular we examine separately the 
rotational and translational channel of relaxation, and 
how their relative contributions change with decreasing 
the temperature. The Letter culminates with a compu- 
tation of the time correlation function as a sum of the 
simple processes with a-priori determined weights , in 
agreement with simulations at different temperatures. 

The model we employ here has particles of mass 
m interacting via the potential 



i ) = U( nj ) + AU(r ij ,e i ,6 j ) 



(1) 



where 
i and 



is the distance between the two particles 
The first term on the RHS of Q is 



the standard isotropic Lennard- Jones potential J7, 
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whereas the anisotropic part of 



the potential is given by 



(2) 



h(x) 
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= (1 -x z f for \x\ < 1 ;h(x) = Ofor|x| > 1 . 

>j) is the included angle between the relative 
= r, — rj and a unit vector Ui (uj) (referred 
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FIG. 1: (Color online). The relaxation function Cn(t) as 
a function of linear time (upper panel), and of logarithmic 
time in the lower panel. The leftmost curve (in red) pertains 
to T = 0.3 and in order to the right, the temperatures are 
T = 0.25, 0.22, 0.20, 018 and 0.17. Note the extreme slowing 
down in this range of temperatures. 



to below as 'spin') which represents the orientation of the 
axis of particle i (j). The function h((9 — 6o)/0 c ) (with 
6q = 126° and 8 C — 53.1°) has a maximum at 6 = Oo, and 
thus 9q is a favored value of Oi. Thus the anisotropic term 
in the potential favors structures of five-fold symmetry. 
The parameter A controls the tendency of five-fold sym- 
metry, and therefore of the frustration against crystalliza- 
tion. The units of mass, length, time and temperature 
are m, a, t = a^J m/e and e/ks, respectively, with ks 
being Boltzmann's constant. According to the numeri- 
cal simulations presented in 0], for A < 0.6 this system 
crystallizes upon reducing the temperature. The ground 
state crystal has an elongated hexagonal structure with 
anti- ferromagnetic ordering of the spins Ui. For A > 0.6 
the system fails to crystallize upon cooling, at least for 



2 



the simulations times reported in 0]. 

The relaxation function: The time auto-correlation 
function of interest was introduced in 0] in terms of the 
spins, in the form 



Cfl(f) = (i/i\o5>*i(*) •«*(<))) 



(3) 



Our own numerical results for this function are shown in 
Fig. [1] As is customary in the field [3| the measured cor- 
relation function was fitted in 0] to a stretched exponen- 
tial form Cij(t) oc exp[— (f/r a )^]. For A = 0.6 the relax- 
ation is of Arrhenius form with a constant value (3 0.95 
for T > T m = 0.46, but (3 was fitted separately for every 
temperature T < T m where it decreases with tempera- 
ture. The relaxation times were fitted to a Vogel-Fulcher 
law t q = To exp[DTo/(T — To)]; below we will show that 
we can reconstruct the correlation function from elemen- 
tary processes with temperature-independent relax- 
ation functions (i.e. without needing to fit anything at 
different temperatures). As stated previously [1, [B|, 0], 
the Vogel-Fulcher fit is deeply misleading in indicating 
a finite temperature where the relaxation time diverge. 
We do not expect any singularity here for any dynamical 
process for any T > 0. 

Statistical Mechanics: To reconstruct the correla- 
tion function from elementary contributions we recall 
how statistical mechanics is constructed for this system 
0. In Fig. [2] we present the three potentials between two 
particles, depending on the orientation of their spins rela- 
tive to the inter-particle vector distance: lowest in energy 
(in blue continuous line) is the case for which both have 
a favored spin orientation; middle, in green dashed line 
(high, in red dash-dotted line) is the potential when one 
(none) of the spins are in a favored orientation. One sees 
that the minima of these potentials occur with signifi- 
cant gaps in their energies, allowing us to now measure 
the average energy of pairs of particles as a function of 
temperature. These averages fall in three distinct ranges, 
such that the range of variation of each energy is much 
smaller than the gaps between the energies, see inset in 
Fig. [5J We denote the three effective energies below 
as 2Eb 2E g and 2E r respectively. Next, Ref [8j defined 
quasi-species in the form of an n-stars, each of which is 
a given particle decorated by the n inter-particle vector 
distances (edges) to its n neighbors, see for example Fig. 
[3J Each such edge is colored according to the spin ori- 
entations. We denote by k the number of red, green 
and blue edges such that n = i + j + k. It turns out that 
in the temperature range of interest (0 < T < 0.5), in an 
overwhelming majority of n-stars (more than 98%) the 
central particle has a spin orientation that is favorable 
with respect to two of its edges, leading to the constraint 
(cf. 0) 



^2(j + 2k)cijk = 4 

ijk 



(4) 
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FIG. 2: (Color online). Potential curves for particle-pairs 
with two spins, one spin or no spin in favored position (blue 
continuous line, green dashed or red dash-dotted line respec- 
tively). Inset: the measured energies of particle pairs, falling 
in three distinct ranges with gaps between them. 




FIG. 3: (Color online). An example of an n-star with n = 
5,i — 2, j — 2 and k = 1. The central particle has a spin 
with favored orientation with respect to edges 1 and 2. Thus 
these edges can be either blue or green, and this central spin 
cannot be favored with respect to any other edge. In the 
interesting range of temperatures we observe 36 n-stars with 
4 < i + j + k < 6. 



which is important for the statistical mechanics of this 
system. The energy of an n-star (referred to as a quasi- 
particle) is computed as 



Eijk = iE r + jE g + kE b 



(5) 



where k < 2. Note that since the energies on the RHS of 
Eq. ([5]) depend on temperature, so does the energy of the 
quasi-particles. Notwithstanding, in the interesting tem- 
perature range the temperature dependence is weak; we 
take the energies of the quasi-particles as T-independcnt. 
The degeneracy gijk of the energy level (number of quasi- 
particles with the same energy) was computed in Q in 
the form 



Qijk 



i + j + k - 
i 



2 J 



+2k-2 A i-k+2 



(0) 
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Finally we can write the partition function of the sys- 
tem: 



Z(T,X(T)) 



ijk 



(7) 



The Lagrange multiplier A is introduced to insure that 
the constraint ((4]) is satisfied. In terms of the partition 
function the mol-fraction of quasi-particles is 



9ijk 



e 0E ijke -\(j+2k) 



(8) 



Z(T,X(T)) 

In Q it was shown that the prediction of this formula 
agrees well with simulations until the system is 'jammed', 
or more correctly, until the simulation time is too short 
to allow the system to equilibrate. Here we will use this 
theory to compute the spin auto-correlation function. As 
noted in we do not need to discuss separately all the 
36 quasi-species, it suffices to bunch groups of to- 
gether into groups with k = 0,1,2. This bunching is 
natural since it stresses the different environments (po- 
tential barriers for unit vector orientation ) seen by the 
quasi-species in each group, helping us to disentangle the 
competing dynamics leading to relaxation. 

Spin decorrealation: we first note that a spin can 
change its angle with respect to an interparticle vector 
distance either due to the spin rotation with the interpar- 
ticle vector distance fixed, or due to translation in which 
the interparticle vector distance changes, cf. [3J Fig. 6. 
We assert that to a good approximation these relaxation 
channels are independent and competing, allowing us to 
write the relaxation function for each group k = 0, 1,2 of 
quasi-particles in a product form: 



t 



t 



T rot (fc,T)' T tI (k,T) 



) ~ /rot(" 



t 



-)x/ tr (- 



t 



'T Tot (k,T) n T (k,T) 
(9) 

where the subscripts 'rot' and 'tr' stand for the rota- 
tional and translational decay channels respectively. We 
reiterate that these functions should be temperature in- 
dependent except through the temperature dependence 
of the typical relaxation times T m t{k,T) and r tr (fc,T). 

Modeling the translation channel is the same for all 
the k groups, since the slowing down of translation in 
this model is dominated by the decrease in concentration 
of the k = quasi-species as the temperature decreases 
0. Indeed, with the concentration of Ck=o(T) decreasing, 
one defines a typical increasing scale £ (T) = 1/ y/ck=o(T) 
whose physical meaning is the typical length of the coop- 
erative event that results in any appreciable translational 
motion. Accordingly the typical relaxation time associ- 
ated with translation grows like 4] 



), 



T tI (k,T)=T° T (k)e 



(10) 



where fi is the characteristic free energy per particle in- 
volved in the cooperative event of length £ and r tr is an 



attempt time (of the order of unity) that may depend 
on k. The relaxation function ftr(t/T tv (k)) is a simple 
exponential, 



f tl (t/T tl (k,T)) 



= e -t/r u {k,T) 



(11) 



and all the non-Arrhenius dependence comes from the 
dependence of £ on T. 

The rotation channel calls for more scrutiny, since the 
k = and k = 1 quasi-species differ significantly from 
the k = 2 quasi-species, as the latter tend to aggregate 
in clusters. We therefore need to deal with them differ- 
ently. Quasi-particles with k = 0, 1 are relatively free 
to rotate, and all that they need to do is to overcome 
the potential barrier for rotation. We thus expect their 
relaxation times to be of simple Arrhenius form, i.e. 



t(*0 



AE Iot (k)/T 



for k = 0, 1 



(12) 



We expect AE YOt (k) to be of the order of a bond energy 
in both cases, which is indeed what is found, cf. Table [U 
The functional form of the relaxation functions turns out 
to be stretched exponentials, i.e. 

U(t/r rot (k,T)) « e -('/WcT))«*> ; for k = Qj 1 

(13) 

The values of the temperature independent parameters 
are shown in Table HI 

For k = 2 we cannot expect such a simple model to 
hold. The reason is that at lower temperatures the k = 2 
quasi-particles aggregate inside clusters whose average 
size increases when the temperature decreases Q. The 
distribution of cluster sizes depends on T, and to repre- 
sent the rotational relaxation function one needs to de- 
compose it into cluster contributions [a, [Z( • Each cluster 
may decay simply with Arrhenius decay time, but the 
ensemble is expected to show a strongly non-Arrhenius 
relaxation time, as shown in 0,0]. For the sake of brevity 
we do not attempt here to derive the rotational relax- 
ation function of the k = 2 quasi-particles but we simply 
measure it to find 

/rot(t/r rot (fc = 2,T)) « exp[-(Vr rot (fe = 2,T)) ' 6 ] , 

(14) 

with the measured value of T mt {k — 2,T) as shown in 
Fig. [4] Note that this relaxation function is again 
temperature-independent (except through the depen- 
dence of r rot (k — 2)) 

We thus posses models for all the elementary relax- 
ation functions in addition to Eq. ()14[1 . The parameter 
used are temperature independent, and are summarized 
in table[l] All that remains is to sum up the contributions 
together with the right weights. 

Summing up the three contributions. Having 
modeled the spin de-correlation formulae due to the two 
channels of decay for the three natural groups of quasi- 
particles, we should be able to predict the total relaxation 
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FIG. 4: (Color online). The logarithm of the relaxation time 
T"rot(fc = 2) as a function of inverse temperature 



TABLE I: Parameters used in the model 





k = 


k = 1 


k = 2 




0.57 


0.45 


0.6 


AE I0 t(k) 


1.27 


1.33 






0.27 


0.32 


0.39 


r rot 


0.24 


2.48 




rg(k) 


2.48 


1.6 


1.07 



function simply by summing up the three contributions, 
with each one weighted by the predicted concentration 
of the appropriate group of quasi-species with a given k. 
In other words, we should plot 

2 

C fl (i) = $> fc (T)F fc (i) , (15) 

fc=0 

and compare the results with the data. This comparison 
is shown in Fig. [5j 

We re-iterate that all the temperature dependence in 
these plots stems from the dependence of the typical de- 
cay time of each elementary process and from the tem- 
perature dependent weights of the concentration of k 




FIG. 5: (Color online). Comparison of the model relaxation 
function with data from the simulations 



groups of quasi-species. The apparent change of j3 in the 
stretched exponential fit to this data Q as a function of 
temperature stems solely from the change of importance 
of the various k groups as a function of the tempera- 
ture, and from the change of importance of rotation vs. 
translation (remember that the translational elementary 
contributions appear here always as pure exponentials). 



The bottom line of this analysis is two- fold; first, we 
demonstrate that it is possible to break apart relaxation 
processes in glass forming system to more elementary 
contributions, understanding much better the origin o 
the complex time dependence of relaxation functions. In 
this point we simply extend what was proposed in @, 0] . 
Second, we stress the danger of straight numerical fits to 
relaxation functions; when these change their functional 
form with temperature, the reason should be sought in 
the existence of a composite process with various dynam- 
ical contributions, each of which may be quite simple. 
Whether the system is 'fragile' or 'strong' in the Angell 
parlance Q may be in the eyes of the beholder, espe- 
cially if the said beholder did not reveal the details of 
the physical phenomenon. 
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